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Abstract: Two complementary analyses of a cyclic negative feedback sys- 
tem with delay are considered in this paper. The first analysis applies the 
work by Sontag, Angeli, Enciso and others regarding monotone control sys- 
tems under negative feedback, and it implies the global attractiveness towards 
an equilibrium for arbitrary delays. The second one concerns the existence 
of a Hopf bifurcation on the delay parameter, and it implies the existence of 
nonconstant periodic solutions for special delay values. A key idea is the use 
of the Schwarzian derivative, and its application for the study of Michaelis- 
Menten nonlinearities. The positive feedback case is also addressed. 

Key Words: delay systems, Schwarzian derivative, Michaelis-Menten func- 
tions, negative feedback, Hopf bifurcation. 

Consider the cyclic nonlinear system 



X n fjn(%l(t 7")) fi n X n , 

for n > 1. Assume that each function ^ is either increasing or decreasing 
and that the system is subject to negative feedback. More formally, let 

fM > 0, Sig'^x) > 0, Si e {1,-1}, i = 1 . . . n, and St ■ . . . ■ 8 n = -1. (2) 
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This system can be considered a generalization of classical models, by 
Goldbeter [JJ and Goodwin |2] , of autoregulated biochemical networks under 
negative feedback. Delay systems with this general structure can also be 
found in the modeling of neural networks, for instance in 011], using Qi{x) = 
a>i tanh(/3jx) as nonlinearities. It should also be noted that different delays 
can be introduced in the nonlinear terms of each equation without loss of 
generality, since all but one of them can be removed with a simple change of 
variables. 

An important special case in biochemical models is that in which those 
functions gi{x) which are not linear have the Michaelis-Menten form 

ctx m d 
f(x) = +c, or fix) = +c, a, b > 0, c > 0, m = 1, 2, . . . . (3) 

w b + x m b + x m ~ 

A recent (though undelayed) model within this framework is that of the 
so-called repressilator, see Elowitz and Leibler jSj. We will give especial 
attention below to this type of nonlinearity. 

The dynamics of the bounded solutions of system (JJJ) under assumptions 
(J2J) is governed by a Poincare-Bendixson result, proved by Mallet-Paret and 
Sell in 1996 Informally speaking, for every initial condition the solution of 
the system approaches either an equilibrium, a periodic orbit, or a homoclinic 
chain of orbits. In particular, any chaotic behavior is ruled out. In the 
positive feedback case 8\ ■ . . . ■ 8 n — 1, system ([!]) is monotone and also falls 
within the framework of Mallet-Paret and Sell. A large number of results are 
known in that case, the most important one perhaps being that the generic 
solution is convergent towards an equilibrium [7||H]. 

The work of Sontag and Angeli [5] can be used to establish a relationship 
between the system (jj) and the one- dimensional discrete system 

u k +i = g(u k ), (4) 

where 

g(u) := —g x o —g 2 o . . . o —g n . (5) 

Namely, if the discrete system (jlj) is globally attractive towards its unique 
equilibrium x , then the original system (JTJ) is globally attractive towards its 
unique equilibrium, for all values of the delay r; see also [TUJ HU C2 PEH IT3]. 
and Hale and Ivanov [To] . 
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A second branch of study for systems analogous to JIJ is the search for 
nonconstant periodic oscillations. This usually involves a different kind of 
assumption, namely that the system (0) is 'ejective' around its unique equilib- 
rium for large enough delay. Such arguments usually require the hypothesis 
|g'(a;o)| > 1, which in particular rules out the global attractiveness of (HJ). See 
Nussbaum [IE] . Hadeler and Tomiuk |17| . Hale and Ivanov [ISJ, and Ivanov 
and Lani-Wayda ^H], among others. 

In the present paper, both approaches are unified to give a more com- 
plete picture of the relationship between system (JTJ) (under assumptions ©) 
and system (0J). A Hopf bifurcation approach is considered to prove that 
|g'(:£:o)| > 1 implies the existence of periodic solutions of (JTJ) for certain val- 
ues of r. Also, an important class of nonlinearities is shown to be such 
that the following conditions are a dichotomy: 

1. The system is globally attractive towards x . 

2. \g'(x )\ > 1. 

The main result of this paper is the following theorem, where Sg denotes the 
Schwarzian derivative of the function g; see for instance |19j . 

Theorem 1 Consider a system under assumptions and let the func- 
tion g(x) defined by be bounded. Suppose that all nonlinear functions gi(x) 
are of Michaelis-Menten form |3J] ; m > 1, or that Sg < 0. Then exactly one 
of the following holds: 

1. System (Q) is globally attractive to a unique equilibrium, and (QJ) is also 
globally attractive to a unique equilibrium, for all values of the delay r. 

2. System O contains nonconstant periodic solutions, and system (QJ) is 
subject to a Hopf bifurcation on the delay parameter r. In particular, 
(QJ) contains nonconstant periodic solutions for some values of r. 

The Hopf bifurcation in the second case is not shown to be supercritical, 
although this seems to be the case from numerical simulations. Bounded 
decreasing functions with negative Schwarzian derivative include — tan _1 (x), 
— tanh(x), and e~ x (on as well as their decreasing compositions. They 
also include the Michaelis-Menten functions above for m > 1, as is shown 
here in Lemma 121 
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Figure 1: Typical solutions of a) system (JTJ) and b) system (J3J), where n — 3, 
9i = 92 = g3 = tan~ 1 (x), = 0.11, /i2 = 2.5, /X3 = 4, and r = 80. c) The 
induced decreasing function and the increasing function g 2 (x) = g(g(x)) 
(see LemmaOJ). Here |g'(xo)| = 1/1-1 < 1- 
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Figure 2: The same system is considered as in Figure ^ except that the value 
of iii has been changed to 0.09. The typical solutions of a) system (JTJ) and 
b) system (JIJ now appear to be limit oscillations and periodic 2-cycles. c) In 
this case |g'(o;o)| = 1/0.9 > 1 and g 2 (x) = x has several solutions. 
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It is important to note that this information is not provided a priori 
by the Poincare-Bendixson theorem itself, which doesn't give conditions for 
the different possible outcomes. Even knowing that the equilibrium of ((TJ) 
is unstable doesn't guarantee the existence of periodic oscillations, since for 
instance homoclinic orbits need to be ruled out (possibly using Morse de- 
composition theory |20j). 

A corresponding theorem will be stated and proved in the positive feed- 
back case 8% • . . . • S n — 1, simplifying a result from Angeli and Sontag [§| in 
the case without delays, and from Enciso and Sontag [21] in the delay case. 

System (@J is one- dimensional and doesn't contain delays, which makes it 
much more tractable than JI}. The assumption of the negative Schwarzian 
is a common simplifying hypothesis in the discrete systems literature; see for 
instance [22] f° r an application to continuous systems. A Hopf bifurcation 
approach has also been proposed in the Poincare-Bendixson context in [25] . 

The direct contributions of the present paper are i) to show that for an im- 
portant class of nonlinear it ies the two alternative cases form a dichotomy; ii) 
to formally establish a relationship between the discrete and the continuous 
system, which has already been conjectured by Smith [21] in the undelayed 
case; iii) to carry out a direct Hopf bifurcation analysis of the linear system 
associated to (JTJ) (which is new to my knowledge), and iv) to illustrate the 
usefulness of the Schwarzian derivative in the context of Michaelis-Menten 
functions. 

In Section n the concept of the Schwarzian derivative is briefly introduced 
and applied to Michaelis-Menten functions. In Section [21 the discrete system 
and its relationship with (0) are described. In Section El the Hopf bifurcation 
argument is developed, and in Section [U the positive feedback case is con- 
sidered. Finally, in Section El the relationship with the general results in [U] 
and [Tl| is shortly discussed, and a conjecture is described from numerical 
simulations. 



An important concept related to the stability of discrete dynamical systems 
is the so-called Schwarzian derivative Sf of a real function /, defined by 



1 Sg and Michaelis-Menten Functions 
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The properties of Sf that will be useful here are summarized in the 
following lemma; see [T^], Section 2B for proofs and details. Intuitively, the 
condition Sf < restricts the form of the function / so that the dynamics 
of Uk+\ = f(uk) is more easily determined. 

Lemma 1 Let f, g be C 3 real functions on a real interval. Then the following 
hold: 

1. If Sf < 0, then f cannot have positive local minima or negative local 
maxima. 

2. S(f o g){x) = Sf(g(x))g'(x) 2 + Sg{x). 

3. Sf < 0, Sg < imply S(f o g) < 0. 

It is now shown that the class of functions with negative Schwarzian 
derivative includes the cooperative Michaelis-Menten functions with m > 1, 
and that S(x/ {b + x)) = 0. 

Lemma 2 Let a, b > 0, c > 0, and m = 1, 2, . . and define 
fix) = h c, gix) = h c. 

77T, 2 — 11 

T/ien Sf(x) = Sg(x) = - 



2 x 2 

Proof. Noting that the Schwarzian derivative doesn't change after multipli- 
cation by or addition of a constant, we can assume that a — 1, c = 0. Using 
the quotient rule we compute 

mx m ~ x mx 2rn ^ 1 m 2 m . 
o + x m (b + x m y x x 

where y = f(x). Similarly we compute 
j ji 

f"(x) = —My - l)(2my - (m - 1)) 



f"'{x) = —^y(l — y) [Qm y + (6m — 6m )y + (m — l)(m — 2)1. 

x A 
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We calculate the Schwarzian derivative 

!/(,)-««) 



[6m 2 y 2 — 6m(m — l)y + (m — l)(m — 2)] 



x 2 



f(x) 2\f'(x) 

i 2 

3 1 

-- — [4m 2 ?/ 2 — 4m(m — l)y + (m — l) 2 ] 
_[( m _ i )(m _ 2 ) - ( m - I) 2 ] = 



x 2LV /v y 2 V ' J 2 x 

To compute Sg(x), it is easy to see that g = b~ x foK, where k(x) = b x l m jx. 
A simple computation shows that Sk = S(l/x) =0, i^O. Therefore 

Sg( X ) = S(foK) = Sf{K(x))K'(x) 2 +SK{x) = -"^Ix 2 1 1 " ~ _ 1 1 



2 ~ x 4 2 x 2 ' 



2 The Discrete System 

Consider a continuous, bounded, decreasing function g : I — > I, where J = E 
ot I = [a, oo), ael. It can be easily seen that there is a unique fixed point 
Xo of g. The study of the discrete system (J3J) becomes straightforward by 
relating its dynamics to that of the system = g 2 (w fc ), since the function 
g 2 (x) = g(g(x)) is bounded and increasing. We state the following lemma for 
convenience; see also Angeli and de Leenheer j2H] for an extended discussion. 

Lemma 3 System is globally attractive if and only if the equation g{g{x)) = 
x has the unique solution xq. 

Proof. All solutions of the system u^+i = g{g{uk)) are monotonic increasing 
or decreasing, and each converges towards some fixed point by boundedness 
and continuity. Furthermore, this system is globally attractive to xq if and 
only if (plj is globally attractive to xq. The conclusion follows immediately. I 

Let I = R or I = [a, oo) and let g : I — > / be differentiable, bounded and 
decreasing. We say that system (JU) is fix point determined if 

Ifl^d)! system (jlj) is globally attractive towards xq. 
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Thus, the global attractiveness of (J3J) is determined by the slope of g{x) at 
its unique fix point. For instance, it was shown in ^2] that the functions 
g(x) = Af(K + x), x > 0, form fix point determined systems for every 
A, K > 0, since for such functions system (jlj) is globally attractive and 
|<7'(xo)| < 1; see also Corollary [TJ 

An example of a (discontinuous) function which is not fix point deter- 
mined is 

{1, x < -0.5, 
0, -0.5 < x < 0.5, (6) 
-1, x>0.5. 

This function has the unique fix point xq = and g'(0) = 0, but there is 
the obvious stable cycle g(l) = —1, g{— 1) = 1. To obtain a proper example 
of a differentiable function which is not fix point determined, it is sufficient 
to smoothen g(x) above with an appropriate convolution operator. 

The reader will have noticed the importance of g being fix point deter- 
mined from the discussion leading to the statement of Theorem ^ Never- 
theless g is only defined in terms of the functions gi, and the composition of 
fix point determined functions is not necessarily fix point determined (nor is 
the composition of merely sigmoidal functions necessarily sigmoidal). This 
is why the Schwarzian derivative becomes useful at this point. 

Lemma 4 Let g : / —>■ I be C 3 , decreasing and bounded, and such that 
Sg < 0. Then g is fix point determined. 

Proof. Consider the increasing function G = g 2 = gog, and note that G'(xq) = 
g'(xo) 2 . If |g'(a;o)| > 1, hence G'(x ) > 1, then by boundedness it must follow 
that G(z) = z for some z > Xq. Therefore (@J) has a nontrivial cycle of period 
2, since g(z) ^ z. 

Conversely, let G'(xo) < 1, and assume that G(z) = z for some z ^ xq. 
Without loss of generality we can assume that G{y) = y, G(z) = z, for some 
y, z such that y < x < z. We use here the fact that SG < 0, by Lemma [T] 
Suppose first that G(x) = x on some interval containing xq. Then in the 
interior of this interval it would hold G"(x) = G"'(x) = and thus SG = 
for these points, a contradiction. It is easy to conclude, using the mean value 
theorem, that there exist constants yi,zi such that y < y± < xo < Z\ < z 
and G\yi) > 1, G'(zi) > 1. Now consider the function G'{x) on the interval 
[yx, Zi]. The results above imply that this function has a minimum ui\ on the 
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interior of this interval, and that therefore G"{w\) = 0, G'"(wi) > 0. Thus 
SG(wi) > 0, a contradiction. I 

Corollary 1 Let I = R or I = [a, oo), and let g : I I be decreasing 
and bounded. If g(x) is the composition of functions each of which either i) 
has negative Schwarzian derivative, or ii) is of Michaelis Menten form for 
m > 1, then g is fix point determined. 

Proof. If g is the composition of functions all of which have negative Schwarzian 
derivative, then this must be true of g as well, and g is fix point determined 
by Lemma 0J The same holds if some of the ^ are of Michaelis- Menten 
form with n > 1, by Lemma El If some but not all of these functions are 
of Michaelis Menten form for m = 1, then still Sg(x) < by the derivation 
formula in Lemma H 

Finally, if all the functions are of the form (ce+ 0x) / (j +5x) , a, f3, 7, 8 > 0, 
then g and g 2 are also of this form. It is then easy to show that the (bounded, 
increasing) function g 2 (x) is concave down on J, and that it has a unique fixed 
point xq which further satisfies g'(xo) 2 = (g 2 )'(xo) < 1. The result follows 
from Lemma El ■ 

The relationship between the nonlinear system and the discrete system 
(J3J) becomes clear in the proof sketch of the following well-studied result. See 
Angeli and Sontag [^j and Enciso, Smith, and Sontag jTTj for an abstract 
formal treatment, as well as Sontag ^3] for a discussion of the embedding 
argument. The use of the lemma by Dancer in this context is new. 

Proposition 1 Consider a system (OJ) under assumption and let g(x) 
be defined by (0). If Op is globally attractive towards Xq, then (OP is globally 
attractive towards a unique equilibrium. 

Sketch of Proof : An elegant result of Dancer |2E1 shows that in an abstract 
monotone system with bounded solutions and a unique equilibrium, all solu- 
tions must converge towards this equilibrium (the result is stated for discrete 
systems in [26 , but a variation for continuous systems is straightforward). 
Consider the extended In- dimensional system 

J- i 
X n 

Zi 



gi{Xi+i) - i = 1. ■ - n 

~- g n {zi{t - T)) - fi n X n , 

gi{zi+i) - fMZi, i = l...n- 
-- gJxAt - t)) - n n z n . 



(7) 
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It is not difficult to see that a trajectory (x\(t) . . . x n (t)) is a solution of (JTJ) if 
and only if (xi(t) . . . x n (t),x\(t) . . . x n (t)) is a solution of (J7J). Moreover, this 
system is now subject to positive feedback, since 5\ • . . . 8 n • S\ ■ . . . ■ S n — 1. 
Thus this system is monotone with respect to a certain partial order; see [S] , 
Chapter 5, and [2Z1- Finally, the equilibria of this system are in bijective 
correspondence with the solutions of g(g(x)) = x. The conclusion follows by 
the result by Dancer and Lemma El 



3 Hopf Bifurcation 

In this section we consider the linearization 

Xi ^{Xi-\-i fXiXij % 1 . . . Tl 1, 
Xn k n X\(t 7") fJj n X n , 

of system (JTJ) around its unique equilibrium point (x~x, . . .x n ). It is easy to 
see that 



(8) 



h = 9i(%i+i), i = l-..n 
kn = g' n (xi). 



(9) 



We will show in the negative feedback case k\ . . . k n < that for \k\ ■ . . . ■ k n \ > 
Hi ■ . . . ■ fi n , a Hopf bifurcation exists on the parameter r. The characteristic 
polynomial associated to the linear system (JHJ) is 

g(z, t) := (z + f i 1 ){z + fi 2 )-...-{z + fi n ) + Ke- TZ , (10) 

where K := — k% • . . . • k n > 0. See Lemma 3 of Hofbauer and So 



Lemma 5 Let g(X, tq) = for A = a + uji, tq > 0, and assume that a > 0. 
Then there exists an open neighborhood U of To, and a differentiable function 
p : U — » C, such that 5f(p(r),r) = on U. If a = 0, then Re p'(r ) > 0. 

Proof. Define f(z) := Yli( z + A^)- The proof of the first statement follows 
by the implicit function theorem for the function g(z, r) at the point (A, To), 
after verifying that dg/dz ^ at that point: 

^(A, r ) = /(A) TT — " ToKe~ XT ° = -K e -^Q(X, r ), 
oz A + Uj 

j J 
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where 



Q(A,r ) ^Vy-j— +r . 

3 J 

Using the fact that > for every j, it is easy to see that Re Q(X, r ) > 
and the proof is complete. 

To prove the second statement, let a = 0. Note that necessarily uj ^ 0, 
since g(z, r) > whenever z > 0. Assume uj > 0, the other case being similar. 
Multiplying on both numerator and denominator by A — fij, it follows that 



Im Q(X, t ) = -uj 2 ]_ 2 < °- 



^ UJ Z + ^, 

3 J 



By the implicit function theorem, 



P'(ro) = "^( A ^o) (|f (A, r ) ) = -uiQ{\T Q y\ 
It follows that Re p'(tq) > as stated. 



Theorem 2 If K > \i\ ■ . . . ■ \i n , then system (GP has a Hopf bifurcation on 
the parameter r. 

Proof. 

We show that there exists r > such that 

i) g(uji, To) = for some uj > 0, 

ii) ^(A, r ) ^ 0, for all A G C with Re A > 0, 

hi) for some uj > 0, it holds that g(ujo,r ) = and that if g(X,T ) = 0, 
A = mujo for integer m then A = ±a;o^- 

Together with Lemma EJ this will directly imply the existence of a Hopf 
bifurcation at the point t = to\ see Theorem 11.1.1 of Hale [29J. 
Let 

S := {t > | g(X, t) = for some A G C such that Re A > 0}. 
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To see that S is nonempty, first note that whenever u > and \ f(u)i) \ = K, 
one can find r > such that e~ MT = —f(ui)/K and so g(ooi, r) = 0. Noting 
that |/(0)| = fix -...-//„< K and \f(ui)\ — > oo as — ► oo, it follows by the 
intermediate value theorem that |/(o;i)| = If for some a;; therefore 5 7^ 0. 

Let To := inf S; it is shown now that r G S. Let 01 > <7 2 > • • • — > To, and 
let Ai, A 2 , • • • be such that Re Aj > and g(\i, <7») = for every z. Let M > 
be such that |/(z)| > 2iT for \z\ > M. Then |e~ AiT | < 1, and therefore 
necessarily |Aj| < M, for every i. There exists thus a subsequence of {Aj} 
which converges towards Ao G C, Re Ao > 0. By continuity g(\o, To) = 0, 
and t G S 1 . 

To complete the proof of i) and ii), it suffices to show that g(\, r ) = 0, 
Re A > imply Re A = 0. But this follows directly from Lemma |3J by the 
minimality of To- 

To see iii), simply recall that g(ui,To) = implies \ui\ < M, and pick 
Uo > so that Uoi is a root with maximal norm. ■ 

Note that this result is proved in the context of Theorem 11.1.1 of 
The existence of periodic solutions for certain values r > r follows, but no 
assertion is made regarding their stability. This may nevertheless be shown 
using the above proof, if the asymptotic stability of the equilibrium of (JTJ is 
established for r = r . 

In the particular case r = 0, it is known j^O] that system (jSJ) is asymp- 
totically stable provided that Kj{pt,\ ■ . . . ■ fi n ) < sec n (7r/n) = l/(cos n (n/n)). 
Therefore necessarily tq > in those cases. 

The following proposition establishes a global stability result for the linear 
system (JSJ). 

Proposition 2 Let K > [L\ • . . . Let r > be the Hopf bifurcation point 
as in Theorem® and let t > 0. Then the linear system (0) is exponentially 
unstable if and only if t > To- 
Proof. It was shown in the proof of Theorem |21 that if r < r (r = r ), then 
g(-,r) could have no root A with Re A > (Re A > 0). Therefore for r < r , 
the exponential stability of (jSJ) is ruled out. 

Let 5" be the set of r > such that system (JHJ) is exponentially unstable. 
It follows from Lemma El and the proof of Theorem |21 that (to,to + e) C 5" 
for some e > 0. Assume by contradiction that S' 7^ (r , 00), and let T\ be the 
infimum of (r , 00) — S'. In particular, it holds that T\ > t + e > r . But 
this is once again a violation of Lemma El I 
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3.1 Proof of Theorem [T] 



Now the proof of the main result is complete. 

Proof. Let xq be the unique fix point of g(x). The first case corresponds 
to the situation in which |g'(:ro)| < 1- Since (@J) is fix point determined by 
Corollary Q and Lemma EJ it holds that (@J) is globally attractive to equi- 
librium. By Proposition ^ system (£Q) is also globally attractive towards a 
unique equilibrium. 

In the case |(/(:ro)| > 1, system (JU) must have a periodic solution since 
it is fix point determined. Evaluating g'(x) using the chain rule yields that 
K > \i\ • . . . • n n , and therefore by Theorem |2l a Hopf bifurcation occurs on 
the parameter r. I 

4 The Positive Feedback Case 

We state and prove the corresponding statement in the positive feedback 
case, which follows from the material in (2] and [21]. Given a system (£Q), 
consider the positive feedback hypotheses 

Hi > 0, Sigi(x) strictly increasing, <5j e {1, —1}, % = 1 . . . n, and Sr. . .-S n = 1. 

(11) 

It is easy to see, by setting the RHS of (JT]) equal to zero, that the func- 
tion x G K. — » <f>(x) = (xi,...x n ) defined by x n := [i~ l g n (x), x n _\ : = 
f J 'n-i9n-i( x n)j ■ ■ ■ i x i '■= lh 9i{ x i) is a bijection between equilibria of (JU) 
and equilibria of (JU). In the following result, 'almost every' is meant in the 
sense of measure, more precisely in the sense of prevalence in For other 
notation used in the proof, refer to [21]. 

Theorem 3 Consider the system (QJ) under Ml)) , and let g be bounded and 
have countable equilibria. Then almost every solution of converges to- 
wards some equilibrium x such that g'{x) < 1. Also, almost every solution of 
(QP converges towards some equilibrium <j>{x) such that g'(x) < 1. 

Proof. The first part of this theorem is straightforward: every solution of ^ 
is monotone increasing or decreasing, and it is bounded since g is bounded. 
Therefore each solution must converge towards an equilibrium. But if x is 
an equilibrium of (J2J) such that g'{x) > 1, then it is repelling and no strictly 
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monotone solution can converge towards x. Since g is strictly increasing and 
therefore injective, given any z e /, k > 1 it holds that g fc (<2) = x implies 
2 = x. Thus the only solution that converges towards x is that with initial 
condition x. 

The second statement follows from Theorem 6 of [21], which is a conse- 
quence of results in [HIE]- It follows from the results in ^U], or from Chapter 5 
of jS] after a change of variables, that (JTJ) is monotone with respect to an or- 
thant cone. The strong monotonicity follows using the strict monotonicity of 
the functions Qi{x) in (fTTj) ; see the proof of Theorem 1.3 in [32]. The bound- 
edness of g implies that one of the functions g>j is bounded, and therefore all 
solutions of the system are also bounded. 

We write system (JTJ as the closed loop of a control monotone system by 
replacing g n (xi(t — r)) on the right hand side by g n (u) and by letting h(x t ) = 
x\(t — r). Theorem 6 of [21] implies in this case that almost all solutions 
converge towards equilibrium points e = (x\, . . .x n ) = <p(x) such that either 
the linearization (JHJ) is not irreducible, g' n (xi) = 0, or else g\x) < 1. But any 
of the first two conditions imply that in the linearization (JHJ) it holds ki = 
for some i, and that therefore 

g'(x) = ki ■ . . . kn/ifMx ■ . . . • fJ, n ) = < 1. 

Thus almost every solution converges to an equilibrium e = 4>(x) such that 
g'(x)<l. I 

5 Future Work 

The framework of Angeli and Sontag [H] and Enciso, Smith and Sontag JT] 
describes quite general dynamical systems as the negative feedback loop of 
controlled monotone systems. Sufficient conditions are then given for the 
system to be globally attractive to equilibrium, even in the presence of de- 
lays or diffusion terms. Theorem can potentially be used to extend these 
results to the case of periodic oscillations, as well as to show that the original 
results are sharp in some sense. It is not the first time that this is suggested. 
For instance, Angeli and Sontag [12] have pointed out that if the associated 
discrete system has a 2-cycle, then large enough delays would create the ap- 
pearance of oscillatory behavior (or pseudooscillations) , which in a biological 
system might be as meaningful as proper periodic oscillations. 
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The analysis of the asymptotic behavior of the system considered in this 
paper is far from complete. If the system falls into the second case of the main 
theorem, simulations suggest that for r > tq the system is in fact globally 
attractive towards a unique nonconstant periodic solution. Work towards 
such a result would most likely include the use of the Poincare-Bendixson 
result, for example by finding a Morse decomposition of the system and 
ruling out the existence of homoclinic orbits. 

Finally, note that the need for the assumption Sg < can be traced back 
to the particular approach used to prove the existence of periodic oscillations 
(Hopf bifurcation), in the following sense: if one could prove the existence of 
periodic oscillations of (£[]) based solely on the existence of a stable periodic 
2-cycle of (j3J), then the proof of the main theorem wouldn't have to require 
that g is fix point determined, and the assumption Sg < could be dropped. 
Indeed, it has been observed in numerical simulations that whenever there is 
a stable 2-cycle of (j3J), then there is also a limit cycle of (JTJ for large enough 
t - - even when Sg -ft 0. This has been numerically observed to be also 
true in more complex noncyclic systems in the framework of [TTJ. Note that 
the proof of the existence of periodic solutions would require to abandon any 
obvious use of Hopf bifurcation or ejective fixed point methods, since it could 
not be required that |g'(xo)| > 1- 
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